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Abstract 



^he duality between QCD perturbative series and power corrections recently conjectured by Narison and Zakharov is analyzed. 
O ,We propose to study correlations between both contributions as diagnostics tool. A very strong correlation between perturbative 
D and non perturbative contributions is observed for several observables at zero and at finite temperature supporting the validity of 
the dual description. 
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1. Introduction 

The disentanglement between perturbative and non per- 
turbative effects in QCD has been a major enterprise in 
the last three decades since the sum rules technique was 
first suggested [112] . While short distance radiative correc- 
tions are perturbatively computable and are characterized 
by a smooth logarithmic dependence on the relevant en- 
ergy scale, non perturbative effects manifest as a stronger 
power-like dependence and in terms of vacuum expecta- 
tion values of local and low dimensional gauge invariant 
operators. Although the underlying quark-gluon dynamics 
should determine the relative strength of perturbative and 
non perturbative contributions unambiguously, up till now 
these condensates have been treated de facto as indepen- 
dent parameters, unrelated to the first few terms of the 
perturbative series. During many years renormalons have 
been viewed as a bridge between perturbative and non per- 
turbative physics (see e.g. Ref. [314] and references therein 
for a review). In a recent paper Narison and Zakharov [5] 
have conjectured a quite different scenario, namely a du- 
ality between condensates and perturbative contributions. 
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This duality concerns the properties of large order pertur- 
bative series which involves an expansion in the strong run- 
ning coupling constant a s (Q 2 ), and it establishes that they 
are dual to non perturbative power corrections, i.e. powers 
of (Aqcd IQ) n ■ At the practical level, this means that if one 
considers a short perturbative series then one should add 
the leading power correction by hand. Only when one uses 
long perturbative series there is no reason to add power 
corrections. The confirmation of this conjecture might pave 
the way for practical approaches where condensates are re- 
quired; the need of allowed and forbidden local conden- 
sates would be justified because of a lack of a complete 
perturbative series to all orders. This is quite timely since 
in many applications at zero temperature and finite tem- 
perature the phenomenological need for power corrections 
involving dimension-2 operators is overwhelming. The zero 
temperature example par excellence is given by the heavy 
qq potential, where the string tension, a dimension-2 ob- 
ject not related to any local gauge invariant operator, ap- 
pears as the unequivocal signal of confinement. There have 
been speculations on the appearance of dimension 2 contri- 
butions to the average plaquette [5] . Recent lattice calcu- 
lations [7] computing the average plaquette in lattice field 
theory up to 20-th order do not yet see the onset of renor- 
malon physics; a mild geometric type perturbative series 
which successfully reconstructs the full result, is observed 
instead. 

While high momenta at zero temperature probes the the- 
ory in the asymptotically free region, we note that a paral- 



lei discussion for finite temperature T 3> Aqcd/27t above 
the deconfinement phase transition might be carried out. 
However, the explicit breaking of Lorentz invariance trig- 
gered by the privileged heat bath reference frame makes the 
theoretical discussion much more involved |8|9j . At finite 
temperature the recent discovery of inverse temperature 
power corrections from relatively old lattice data becomes 
evident from plots in 1/T 2 and has been quite impressive 
and rather unexpected. They can effectively be explained 
by a dimension two condensate, starting from the Polyakov 
loop |10|llj , the heavy quark-antiquark free energy in [12] 
as well as the trace anomaly and QCD equation of state 
in [13114115] (see also Ref. [16]). Of course, it would be quite 
interesting to determine whether these thermal power cor- 
rections are dual, in the sense of Narison and Zakharov, to 
a long perturbative series. The present paper addresses this 
important issue. 

While the duality conjecture might eventually be tested 
more convincingly in the future, we suggest an alternative 
approach where some quantitative insight may also be gath- 
ered both at zero and finite temperature from confronting 
current lattice and perturbative results. Basically, the idea 
is quite simple. Given that the only scale entering the per- 
turbative series to any order is Aqcd, in a fit to lattice 
data containing both perturbation theory and condensates 
we should observe 1) smaller contributions from condensate 
at increasing perturbative orders and 2) a rather strong 
statistical correlation between Aqcd and the dimcnsionful 
condensates. 

In the present paper we pursue this idea to check with 
the accessible theoretical information and lattice data the 
validity of the duality conjecture. In Sec. [2] we deal first 
with the more familiar quark-antiquark potential at zero 
temperature, where we indeed observe, within uncertain- 
ties, the expected correlations. This not only supports the 
perturbative-power duality but also qualifies the correla- 
tion method as a handy tool to study the duality elsewhere. 
We are interested to do so at finite temperature above the 
deconfinement phase transition for the Polyakov loop in 
Section [3] and the trace anomaly in Section [4] As a useful 
guideline we use the finite temperature model where non 
perturbative thermal power corrections are driven by a di- 
mension 2 gluon condensate in the dimensionally reduced 
theory [10111112113114115] . While the model might be im- 
proved, taken at face value it provides a unified and coher- 
ent description of gluodynamics lattice data using the same 
values of condensates within estimated errors. It therefore 
provides an ideal playground to search for possible correla- 
tions in the sense of the above mentioned duality. 

2. The quark-antiquark potential 

Our points are best exemplified with the (zero temper- 
ature) heavy quark-antiquark potential within a perturba- 
tive expansion for which many efforts have been devoted. 
Nowadays the perturbative series are partially known up 



to order a* log 2 a s in the Weyl or temporal gauge, Aq = 
(see Ref. [17] and references thereirP~h . In either case it is 
known that this perturbative computation can only repro- 
duce lattice data at small separations, and fails for r > 
0.25 fm. The Cornell potential is a phenomenological ver- 
sion of this potential and it gives a good overall description 
of the lattice data [5D] for all separations. It reads 

V,?(r) = -~ + or, (1) 

where a s is the QCD coupling constant, which is consid- 
ered as a constant, i.e. no running, and a ~ 4.64 fm -2 is 
the string tension term. The Coulomb term corresponds to 
the leading order in perturbation theory. The linear term 
follows from quarkonium phenomenology, and it is widely 
accepted that it cannot follow from a perturbative compu- 
tation. This leads to the conclusion that the perturbative 
potential is not complete, and should be extended with a 
linear term put by hand, i.e. 

1 N 

Vgg(r) = -^2a„a^(r)+a N r. (2) 

Because the accepted separation between perturbative and 
non perturbative contributions, one is tempted to identify 
the parameters o~n and a. In this section it will be shown 
that this identification might not be correct, in line with 
Ref. [5], and in fact there is a mixing between power-like 
corrections and the perturbative series. This means that 
eventually limjv^oo erjv = 0. In order to provide further 
convincing evidence that this might happen, we will analyze 
lattice data for the heavy qq potential from Ref. [5D] using 
Eqs. (D and ©. 

Two kinds of fits are considered. In the first ones, a dis- 
tance interval is sought where the perturbative expansion 
works well. As expected, this corresponds to taking a suf- 
ficiently small distance region, namely, 0.085 fm < r < 
0.170 fm. This regime follows from the requirement that the 
fitted (tjv is compatible with zero within errors. The second 
type of fits use all distances of lattice data, 0.085 fm < r < 
0.830 fm. In both cases an additive constant is allowed in the 
potential, chosen so that the fit reproduces exactly a point 
of lattice data at an intermediate distance r = 0.25 fm. It 
has been checked that the conclusions are unchanged when 
this parameter is also included in the fit. 

For the small distance regime, the perturbative series 
to 0{al) describes very well by itself the lattice data 
(x 2 /dof = 0.04 <C 1) and the string tension cat turns out 
to be compatible with zero. The fits including all lattice 
data are summarized in Table [T] The value of ctjv tends to 
decrease (within errors) as higher orders in perturbation 
theory are included, and this is what would be expected 
from the scheme of duality between power corrections and 
perturbation theory. 



1 See also |18I19| for a recent complete three loop calculations in 
Feynman gauge. 



Order 


S 


cr[fm- 2 ] 


r(c, <5) 


X 2 /dof 


Tree Level 




4.99(11) 


0.991 


0.98 


1-loop 


3.66(7) 


4.25(6) 


0.974 


0.20 


2-loop 


4.54(10) 


4.12(6) 


0.978 


0.79 


N 3 LO |18|19| 


4.31(13) 


4.07(6) 


0.980 


0.93 


N 3 LL Wf] 


4.19(14) 


3.85(7) 


0.984 


0.76 



Table 1 

Fit, using Eq. J2j, of heavy qq potential lattice data from Ref. I20| . 
The fit is performed in the interval 0.085 fm < r < 0.332 fm for the 
tree level, and in the full interval of lattice data, 0.085 fm < r < 
0.830 fm, for the rest. The intervals are chosen so that x 2 /dof < 1. 
At three level, the correlation coefficient refers to — r(<r, a s ). 
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Fig. 1. Correlation ellipses corresponding to A\ 2 = 2.3,4.6, between 
the perturbative parameter <5 = log(/x / Aqce>) an <i the non pertur- 
bative a for the heavy qq potential. The perturbative formula has 
been considered at 3 loops order (N 3 LL) [T7], i.e. N = 4 in Eq. 10. 
The fit uses data in the interval 0.085 fm < r < 0.170 fm. 
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Fig. 2. Same as Fig. [T] but for the last line in Table [T] The fit uses 
data in the interval 0.085 f m < r < 0.830 fm 

The most relevant feature uncovered by this analysis is 
the strong correlation found between the perturbative se- 
ries and the string tension term, not only in the small dis- 
tance regime, but also in the entire distance regime. Also 



the correlation is larger at higher orders[f] Figs. [T] and [2] 
show the correlation ellipses when terms up to N 3 LL are 
included in the potential. This strong correlation confirms 
the dual description of QCD proposed by Narison and Za- 
kharov [5] . In what follows we aim to apply the same cor- 
relation method to determine whether or not this duality 
takes place also at finite temperature above the deconfinc- 
ment phase transition. 

3. The Polyakov loop 

The vacuum expectation value of the Polyakov loop in a 
gauge in which A$ is time independent reads 



L(T) = (— tv c e^ A ^l T 



and an expansion of the exponential gives [lQj : 



log L(T) 



9 2 (Al a ) 
AN r T 2 



0(g 5 )- 



(3) 



(4) 



It has been shown in a series of works [10111112113114115] 
that power corrections provide the bulk of observables at 
finite temperature in the non perturbative regime of the 
deconfined phase of QCD, i.e. in the regime T c < T < 6T C . 
Our considerations were first proposed to describe the lat- 
tice data for the renormalized Polyakov loop in this regime, 
and follows from the introduction of a tachyonic gluon mass 
at short distances in the gluon propagator [10] . This is the 
analog of the zero temperature modification proposed in 
Ref. [21] . Moreover, as shown in [12], a common value for 
this mass reproduces both the string tension and Polyakov 
loop data. These considerations imply that the perturba- 
tive value of (Aq a ) should be augmented with a non per- 
turbative term directly related to the tachyonic mass [10] : 



(5) 



Up to radiative corrections, the perturbative part (Aq a ) p 
is proportional to T 2 whereas (A 2 a ) NP is temperature in- 
dependent. Thus, the total Polyakov loop can be separated 
into perturbative and non perturbative contributions in 
(4J, and reads (N c = 3) 



logL(T) = logL P (T) 



9 2 (Aj a ) 
12T 2 



NP 



(G) 



The presently available perturbative calculations have been 
carried out to order c/ 4 [22] (recently corrected in |23j)f^1 



2 A fit in the interval 0.085 fm < r < 0.830 fm using the tree level 
leads to r(a,a s ) = —0.960. 



For gluodynamics with N c 



logL P (T) = g + |l(log 9+4 



3 this gives 

r 



(7) 



(where the subindex P stands for perturbative). Since the j3 function 
starts at order g 3 (fJ.), changes in /i affect 0(g 5 ). 
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Fig. 3. Correlation ellipses for the Polyakov loop corresponding to 
Ax 2 = 2.3,4.6 between the perturbative parameter S and the non 
perturbative parameter b^. The perturbative Polyakov loop is in- 
cluded to order g 4 . 

Clearly, Eq. ^ , with a perturbative series plus a dimension 
2 power-like term, resembles Eq. 

The lattice data from Ref. [24] for the renormalized 
Polyakov loop with iV c = 3 can be fitted using Eq. ©, 
including the non perturbative term and different orders 
in the perturbative series. The scale in g([i) is taken as 
fi = e s 2nT. Results are shown in Table [5] Within uncer- 
tainties it can be seen that the coefficient of the power-like 
term — g 2 {A\ a } NP /6T 2 does not change much while 
the correlations are quite strong. The corresponding cor- 
relation ellipses arc displayed in Fig. [3J 



Order 


L P (6T C ) 


<5 


b L 


r(b L ,S) 


X 2 /dof 


Lp = const 


1.121(8) 




1.72(5) 


-0.472 


0.45 


0(a 3 / 2 ) 


1.125(11) 


-0.72(20) 


2.23(16) 


-0.957 


1.22 


0( Q 2 ) 


1.123(9) 


-0.06(18) 


2.15(11) 


-0.901 


1.44 



Table 2 

Fit of lattice data for the renormalized Polyakov loop from Ref. 1241 , 
for TV 3 X N T = 32 3 X 8, using Eq. J6} in the interval 1.03 < T/T c < 6. 
The fit is made using the perturbative series Lp(T) up to a given 
order, c.f. Eq. J7jl, in addition to the non perturbative term of Eq. l(6}. 
bij stands for g 2 {Ag a ) NP /6T 2 . In the first line Lp is taken to be a 
constant and the correlation coefficient there refers to — r(fe^,Lp). 

The possible contribution of a dimension four condensate 
in L(T) can also be considered. To this end, we carry out 
a fit with the formula: 



lo g L(T) = logL P (T) 



bL 
2 





f-f 











(8) 



The results are shown in Table [31 It is noteworthy that the 
correlations ripL, 8) tend to increase with the perturbative 
order. Although the fit does not yield a clear signal for cl, 
the correlation r(pL, cl) is remarkably high. This is in line 
with some appreciation in [5] about the eventual extension 
of the duality also to dimension four condensates. 



Order 


L P (6T C ) 


<5 


b L 


CL 


r(b L ,S) 


r(b L ,c L ) 


x 2 

dot 


constLp 


1.118(12) 




1.61(27) 


0.13(34) 


-0.605 


-0.974 


0.41 


0(a 3 / 2 ) 


1.128(16) 


-0.78(28) 


2.39(46) 


-0.14(33) 


-0.845 


-0.871 


1.55 


0(a 2 ) 


1.128(15) 


-0.15(27) 


2.36(44) 


-0.21(37) 


-0.828 


-0.940 


1.61 



Table 3 

Same as table E] but including a term ~ c L /T 4 , c.f. Eq. (§}. The fit 
is made in the interval 1.03 < T/T c < 6. 

4. Trace anomaly and equation of state 



Power-like 1 /T 2 terms have also been found in the equa- 
tion of state of gluodynamics and QCD [16125(15] . The nat- 
ural observable to display such effects is the trace anomaly, 
or interaction measure: 



3P = T*4f£ 



dT V T 4 



1(g) 
2.9 



<(G° ) 2 }- 



(9) 



Unlike the pressure, the trace anomaly gets no contribution 
from the ideal gas part and this is the primary quantity 
used in lattice to obtain the pressure [26] . Perturbatively, 
the quantity ((G^„) 2 } is proportional to T 4 up to radia- 
tive corrections. Within the same model used to analyze 
the Polyakov loop, it has a further power-like contribution 
proportional to T 2 from the dimension 2 condensate [15] : 



<(G" ) 



a \2\NP 



-6™ 2 D (A 2 . a ) 



NP 



(10) 



where mu is the Debye mass. The trace anomaly can then 
be expressed in the following form 



(6 - 3P) (e - 3P) 



pert 



rp4 



(11) 



where 6a = —Zgj3(g)(A^ Q ) NP /r 2 for three colors and no 
quarks. This pattern is similar to the one in Eq. ^ for the 
qq potential and in Eq. ([6]) for the Polyakov loop. 

One can make the same analysis that was performed 
in previous sections, and fit lattice data for the trace 
anomaly from Ref. [26] using Eq. ([TT]) . The weak cou- 
pling expansion for the free energy is known up to order 
fl£log(sr s ) [27128129130131132133134] . A renormalization 
group invariant (RGI) resummation, to be used below, is 
presented in [15] . The perturbative series is poorly conver- 
gent for the lattice QCD available temperatures T < 5T C , 
and even at much higher temperatures. There have also 
been numerous attempts to resum perturbation theory in 
order to get a better convergence of the result, one of the 
most developed techniques being the hard thermal loop 
(HTL) perturbation theory. The free energy of the gluon 
plasma has been computed recently up to three-loop order 
in HTL (see [3S] and references therein). 

In the fit of the lattice data for the trace anomaly we con- 
sider both resummations (RGI and HTL), which enter as 
the term (e— 3P) pert in Eq. (fTTj) . All the fits were performed 
in a regime in which x 2 /dof < 1, so that reliable errors 
could be extracted. In particular, because PT is expected 
to work better as the temperature increases, it is sufficient 
to change the lowest temperature value of the interval. As 



A 



fitting parameters we take 6 a and the parameter S denned 
by /i = e s 2nT. Because the running coupling depends on 
^/Aqcd, a change in 5 is related to a change in Aqcd, and 
the correlation between 6a and Aqcd is measured by the 
quantity t(6a, 5). 



Order 


Ap Cr t/T 4 T=4.5T C 


<5 


&A 


r{b A ,S) 


X 2 /dof 


a 3 =const 
0{a 2 ) 

0(a 5 / 2 ) 
0{a s ) 

0(a 7 / 2 ) 
0(a 4 ) 


-0.02(4) 

04+ ' 04 
u - uy: -0.02 

-0 04+ - 02 
-0 07+ - 03 

u - u ' -0.04 

06+ 05 
U - UD -0.02 

-0.006 ± 0.043 


2.3 ± 1.9 
2.3 ± 1.6 
2.3 ±0.7 
2.3 ±0.9 
0.45 ±2.4 


3.46(13) 
3.29(13) 
3.57(18) 
3.73(19) 
2.74(58) 
3.35(49) 


-0.730 
0.722 
-0.867 
-0.880 
0.891 
0.983 


0.35 
0.86 
0.42 
0.78 
0.99 
0.37 



Tabic 4 



Fit to trace anomaly lattice data from Ref. [26], N% x N T = 32 3 X 8, 
using Eq. lilt . The temperature interval 1.13 < T/T c < 4.54 is used 
for all lines except O(o 7/2 ), for which 2.0 < T/T c < 4.54 is taken. 
In the first line no running in a 3 is applied. In the other lines the 
renormalization group invariant a(T) = 47r/(22 log(2-7rT/AQCD)) is 
used 15 . At 0(a 4 ) the value A 6 = 20.0 is adopted [15]. In the first 
line, the correlation coefficient refers actually to — r(b A ,a s ). 

In the RGI resummation, at order 0(a A ) an undeter- 
mined parameter, Aq, appears due to infrared divergences. 
Setting S = 0, we fit the lattice data using 6a and Aq as 
free parameters at order 0(a A ). In this case, the best fit in 
the regime 1.13 < T/T c < 4.54 gives [15]: 

6 A = 3.18(74), A 6 = 20.0 ± 10.5, r(6 A , A 6 ) = 0.992 , (12) 

with x 2 /dof = 0.40. Next we fit 6 A and 5 from 0(a 2 ) to 
0(a A ). For the highest order the central value of A 6 pre- 
viously obtained is used. The results are shown in Table 2] 
and Fig. 2J The fit favors vanishing small perturbative con- 
tributions, producing large values of e s (> 10 3 ). To prevent 
this, from 0(a 2 ) to C(a 7 / 2 ) an upper bound e s < 10 has 
been set. Although the results are not fully conclusive, in 
general, the correlation between perturbative and non per- 
turbative terms becomes larger when higher orders in PT 
are included, and also, the value of the non perturbative 
coefficient 6a tends to be smaller for the higher orders. The 
correlation ellipses for 0(a 4 ) are displayed in Fig.QJ 

We have also performed the analyses of this section using 
lattice data from Refs. [35] and [37]. The values of the pa- 
rameters agree within estimated errors with those quoted 
here from Ref. [26] . As a rule these alternative lattice data 
lead to larger errors. For instance, the fit using the lat- 
tice data from Ref. [36] in the regime 1.14 < T/T c < 3.6 
yields 6 A = 3.57(54), 3.25(46), 3.72(90), 3.73(80), 3.09(45), 
2.7(1.6), for orders from LO to 0(a 4 ) respectively. The de- 
crease of 6a as the PT order in increased is more evident 
with these data, although they are affected by larger er- 
rors. The correlation increases in this case from r(b/\ , ct s ) = 
-0.741 at LO to r(b A ,5) = 0.975 at 0(a 4 ). 

In order to extract the possible contribution of a dimen- 
sion four condensate, we have also considered a fit using 
the formula: 
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Fig. 4. Correlation ellipses corresponding to A\ 2 = 2.3, 4.6 (dof = 8) 
between the perturbative parameter 5 and the non perturbative 
parameter b A for the trace anomaly (see main text). The perturbative 
series (e— 3P) pC rt is considered to 0(» 4 ) in the RGI resummation of 
Ref. [T5| . The fit uses lattice data in the interval 1 . 13T r <T< 4.54T C . 

(i^ = (i^ +64 (py +C4 (^\ (13) 

The results are shown in Table [5] The parameter ca is 
compatible with zero, but a high correlation between both 
non perturbative terms 6a and ca is displayed, in line with 
results of Sec.[3J The correlation t(6a, 8) tends to increase 
with the perturbative order, as in Table [4] The correlation 
ellipsoid corresponding to a joint fit of A&, 6a and ca is 
displayed in Fig. [5] 



Order 


8 


b A 


CA 


r(6 A ,<5) 


7-(&a,ca) 


X 2 /dof 


a s =const 




3.7(7) 


-0.3(9) 


-0.850 


-0.969 


0.31 


0(a 2 ) 


2.3 ±3.5 


3.09(51) 


0.33 ±0.67 


0.747 


-0.938 


0.80 


0(os>/ 2 ) 


1.6 ±3.1 


3.9 ± 1.2 


-0.5 ± 1.1 


-0.950 


-0.961 


0.30 


0(a 3 ) 


2.3 ±2.1 


4.1 ± 1.1 


-0.6 ± 1.0 


-0.940 


-0.956 


0.31 


0(a 7 / 2 ) 


2.3 ±1.9 


2.6(9) 


1.1 ± 1.1 


0.883 


-0.926 


0.76 


0(a 4 ) 


-0.6 ± 1.2 


2.8 ± 2.0 


-1.2 ±5.5 


0.983 


0.960 


0.29 



Tabic 5 



Same as Table [5] but adding a term c A (T c /T) 4 , c.f. Eq. (T5J. The 
interval 1.13 < T/T c < 4.54 has been used, except at 0{a 7 / 2 ) 
where 1.24 < T/T c < 4.54 has been used. At 0(« 4 ) the value A 6 = 
20.0 ± 10.5 has been adopted, cf. Eq. ([12}. 

The analogous analysis can be carried out considering 
for (e — 3P) pcr t the HTL perturbation theory result at 1- 
loop [38] , 2-loops [39] and 3- loops [35] . The results are pre- 
sented in Tables [5] and [7] 

A large correlation between 6a and <5 (or Aqcd) is found 
even at 1-loop order. In addition, the effect of a smaller 
contribution from 6a at increasing perturbative orders is 
rather clear in the HTL scheme. Once again, the correlation 
between dimension 2 and dimension 4 condensates turns 



r. 



Order A H tl/T 4 |t=4.5T c S b A r(b A ,5) x 2 /dof 

1- loop -O.Q3l°;^ -0.03 ± 0.69 3.69(40) 0.975 0.36 

2- loops — 0.004+g;°o4 -0.42 ± 0.48 3.57(28) 0.949 0.41 

3- loops 0.08(5) 2.3 ±9.7 2.3 ± 1.5 0.977 0.67 
Tabic 6 

Fit of trace anomaly lattice data from Ref. [26], N% X N T = 32 3 X 8, 
using Eq. mil l, in which (e — 3P) pcr t is identified with the HTL 
result [3§l39 -35 . The interval 1.13 < T/T c < 4.54 is used for 1-loop 
and 2-loop orders, and 2.29 < T/T c < 4.54 at 3-loops. 



Order 


<5 


bA 


CA 


r(b A ,5) 


r(&A,CA) 


X 2 /dof 


1-loop 


0.7 ± 1.8 


4.3 ± 1.0 


-0.6 ± 1.1 


0.932 


-0.976 


0.29 


2-loops 


-0.47(80) 


3.58(40) 


-0.1 ± 1.0 


0.411 


-0.145 


0.47 


3-loops 


2.3 ±9.3 


2.37 ± 1.6 


1.3 ± 1.3 


0.964 


-0.898 


0.92 



Table 7 

Same as Table |6] but adding a term c A (T c /T 4 ), c.f. Eq. (15) . The 
interval 1.13T C < T < 4.54T C is used at 1- and 2-loop orders, and 
1.24T C < T < 4.54T C at 3-loops. 



out to be strong0 

5. Conclusions 

We have proposed to study statistical correlations be- 
tween perturbative series and power corrections when ana- 
lyzing lattice QCD results as a way to probe quantitatively 
the duality proposed by Narison and Zakharov in Ref. [5] . 
We have found that the effects of this duality start feeling at 
leading order in perturbation theory, because even at this 

4 The relative small correlation at 2-loops in Table [7] is confirmed 
when other lattice data are used I36I37| . but this result seems to 
be anomalous in view of Table \E\ Nevertheless at this order a high 
correlation between ca and 5, r(cA, S) = 0.840, is found. 




Fig. 5. Correlation ellipsoid corresponding to A% 2 = 3.53 (dof = 7 
and confidence level 68%) between the perturbative parameter Ag 
and the non perturbative parameters 6a an d ca for the trace 
anomaly (see main text). The central values A$ = 29.9, 6a = 4.15 
and ca = -0.53 follow from a fit in the interval 1.13T C < T < 4.54T C . 



order the correlations are very strong. We have performed 
an analysis for the heavy qq potential at zero tempera- 
ture, and for several observables in the deconfined regime 
of thermal QCD, for which power corrections were derived 
in previous works. As a byproduct, we have addressed the 
important question on the finding of inverse power correc- 
tions at finite temperature above the deconfincmcnt phase 
transition. Our analysis dissolves the apparent contradic- 
tion between the real existence of thermal 1 /T 2 power cor- 
rections in the lattice results and the persistent failure of 
perturbation theory to a given finite order to reproduce 
them; these two extremely disjoint scenarios are actually 
complementary and strongly interrelated. 

The present correlation study can be extended to other 
cases where the condensate-perturbative duality might be 
expected. A very extreme situation corresponds to just use 
power corrections and no perturbation theory at all. This 
could be considered the starting point of the analysis, where 
some "non perturbative" physics would be expected. An- 
other extreme situation requires a complete knowledge of 
perturbation theory to all orders, a certainly unrealistic sit- 
uation. For the situation in-between note that to observe a 
decreasing condensate for increasing orders in perturbation 
theory it is not at all trivial as it depends on the behavior 
of the perturbative expansion. On the other hand, show- 
ing that the residual condensate actually vanishes when all 
terms in perturbation theory are taken into account seems 
to us as difficult as solving QCD exactly. 

We acknowledge useful correspondence with the authors 
of Ref. Q3 . 
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